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' Abstract 

, By a similar idea for the construction of Milnor's gamma functions, we introduce "higher 

■ depth determinants" of the Laplacian on a compact Riemann surface of genus greater than one. 

C*") , We prove that, as a generalization of the determinant expression of the Selberg zeta function, 

this higher depth determinant can be expressed as a product of multiple gamma functions and 
, what we call a Milnor-Selberg zeta function. It is shown that the Milnor-Selberg zeta function 

admits an analytic continuation, a functional equation and, remarkably, has an Euler product. 
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1 Introduction 

In 1983, Milnor ([M]) introduced a family of functions {"f r (z)}r>i what he thought as a kind of 
"higher depth gamma functions". These are defined as partial derivatives of the Hurwitz zeta 
function ^(w,z) := X]nLo( n + z )~ w a ^ non-positive integer points with respect to the variable w. 
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We call J r (z) a Milnor gamma function of depth r and will denote it by T r (z) (see [K0W2] for 
several analytic properties of T r (z)). The purpose of Milnor 's study about such functions is to 
construct functions satisfying a modified version of the Kubert identity f{x) = m s ~ l YlT=o /(^r) 
([Kub]), which plays an important role in the study of Iwasawa theory. 

The aim of the present paper is, as an analogue of the study of the Milnor gamma functions, to 
study "higher depth determinants" of the Laplacians on compact Riemann surfaces of genus g > 2 
with negative constant curvature. Let us give the definition of the higher depth determinant in 
more general situations. Let A be a linear operator on some space. We assume that A has only 
discrete spectrum with eigenvalues = Ao < Ai < A2 <■■■—>■ 00. Define a spectral zeta function 
(a(u>, z) of Hurwitz's type by 

00 

( A (w,z) :=X> i + z)- w . 

j=0 

We further assume that the series converges absolutely and uniformly for z on any compact set in 
some right half w-plane, and can be continued meromorphically to a region containing w = 1 — r, 
for r G N. Moreover, we assume that it is holomorphic at w = 1 — r. In such a situation, we define 
a higher depth determinant of A of depth r by 

( d 

Det r (^ + z) := exp (a(w,z) 

\ OW w =l-rJ 

When r = 1, this gives the usual (zeta-regularized) determinant of A. Notice that if A is a finite- 

A r_1 

rank operator and Ai, . . . , Ajv are their eigenvalues, then we have Det r (^4) = rij=i \ 3 > whence 
Det r (,4 © B) = (Det r yl) • (Det r i?) if both A and B are finite-rank. 

To state our main results, let us recall the case of the determinant of the Laplacian, that is, 
the case r = 1. Let HI be the complex upper half plane with the Poincare metric and V a discrete, 
co-compact torsion- free subgroup of ST2(]R). Then, r\H becomes a compact Riemann surface of 
genus g > 2. Let Ar = —y 2 {-j^i + Jp) be the Laplacian on r\M, Xj the jth eigenvalue of Ap and 
put Spec (A r ) := {Xj \ j G N }. We write Xj = r| + \ where rj £ iR>o if < Aj < \ and rj > 
otherwise. Moreover, let := i ± irj. Notice that o£ G [0, 11 with af < if < A, < \ and 
Re {a - ) = \ otherwise. It is shown that the series C/\ T {w,z) converges absolutely for Re (w) > 1, 
admits a meromorphic continuation to the whole plane C and is in particular holomorphic at w = 
(see, e.g., [DHP, S, Vo]). Moreover, the determinant det (Ap — s(l — s)) := Deti(Ar — s(l — s)) 
can be calculated as 

(1.1) det (A r - s(l - s))= G T {s)Z T (s) = 0(s) 9 " 1 Z r (s). 

Here, </>(s) is a meromorphic function defined by 

0( s ) : = e -2(-|) 2 -4(,-i)C'(0)+4C'(-l) r(s) -2 G(s) -4 
= e -2(-i) 2 ri ( a )-2r 2 ( S )4 



with C( s )> L(s), G{s) and r n (s) being the Riemann zeta function, the classical gamma function, the 
Barnes G-function (= a double gamma function ([Bl])) and the Barnes multiple gamma function 
([B2]), respectively. The other factor Zy{s) in (1.1) is the Selberg zeta function defined by the 
Euler product 

00 

Zr(s):= J] H(l-N(P)- s - n ) (Re( S )>l). 

PePrim (r) n=0 
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Here, Prim (T) is the set of all primitive hyperbolic conjugacy classes of T and N(P) is the square 
of the larger eigenvalue of P G Prim(r). It is known that Z-p(s) can be continued analytically to 
the whole plane C and has the functional equation 



(1.2) Z r (l -s)= (Si( fl )^ 2 (s)-*j Z T (s), 

where S n (s) is the normalized multiple sine function defined by (3.31) ([KK]). Notice that if we 
define the complete Selberg zeta function Hr(s) by 

Zr(s) := (^(sfTiis + lfy'Zris), 
then, the functional equation (1.2) is equivalent to 

(1.3) S r (l-*) = S r («). 

Moreover, it is known that Z-p(s) has zeros at s = 1, 0, —k for k G N with multiplicity 1, 2g — 1, 
2{g — l)(2k + 1), respectively (the latest are called the trivial zeros because they also come from 
the gamma factor as the Riemann zeta function) and at s = at for j G N (these are called the 
non-trivial zeros). In particular, Zy(s) satisfies an analogue of the Riemann hypothesis, i.e., all 
imaginary zeros of Zr(s) are located on the line Re (s) = ^. See [D] for arithmetic trial of a 
determinant expression for the Riemann zeta function. 

In this paper, as generalizations of the case of r = 1, we study the function 

Dr,r(») := Det r (A r - s(l - s)):= exp(-^-C Ar (w, -s(l - s)) _^ ) 

for an arbitrary r G N. Here, for j corresponding to the zero and the exceptional eigenvalues A,- 
(that is, < Xj < |), the summand (Xj — s(l — s))~ w of the spectral zeta function (a t (w, — s(1 — s)) 
is defined by (A,- — s(l — s))~ w := exp(— iolog^(Aj — s(l — s))) where log^^ is one of the branch 
of the logarithm defined by (2.4) in Section 2. For the other j, as usual, we employ the principal 
branch of the logarithm log which takes values in IR + i(—ir, ir]. We will see in Section 2 that the 
defining domain of Drr( s ) is a certain region with two connected components which are convex 
and symmetric with respect to the line Re (s) = ^. 
Define the region fir by 

(1.4) fi r :=C\(( II [a|, at - ioo) U [aj,aj + ioo) 



\(( U i 

u ( U [«/' \ + icc ) u K 7 4 " ioo 0) 



Here, for a, /? G C, the semi-closed line from a to are denoted by [a, (3) and so on. Notice that 
VLy is connected if and only if | Spec (Ar) (see Figure 1). The following theorem describes main 
results of this paper. 

Theorem 1.1. Let r > 2. 

(i) The function Dr, r (s) can be extended holomorphically to the region Vi-p and satisfies Dr 5 r(l — 

s) = D r , r (s). 

(ii) It can be written as 

(1.5) D r ,(s) = G r , r (s)Z r , r (s) = M^ZrA 8 )- 
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Figure 2: (p r (s) is holomorphic in C\ ((— oo, — 1]U Figure 3: Zr r (s) is holomorphic in J7r\(— oo, —I] 
[0,-ioo)). (i^Spec(A r )). 
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Here, 4> r {s) is a holomorphic function in C \ ((— oo, — 1] U [0, — ioo)) expressed as (3.16) by the 
Milnor gamma functions and, further, can be written as a product of multiple gamma functions. 
The other factor Zr t r(s), which we call a Milnor-Selberg zeta function of depth r, is a holomorphic 
function in fir \ ( — oo, —1] having a functional equation 

2r -l 
(1.6) Z r>r (l - s) = ([] S»- a ^ s -i)) Z r>r (s) 

j'=i 

for s £ !!r\ ((— oo, — 1] U [2,+oo)) where a r j(t) is the polynomial defined in (3.27). Moreover, 
Zr,r(s) has an Euler product expression in a sense that it can be written as a product and quotient 
of lL poly-Selberg zeta functions" Z r (s) defined by an Euler product for Re (s) > 1. 

Similarly to the result of r = 1, we may call Gr,r( s ) (or <f> r ( s )) a gamma factor. Note that 
4> r (s) and Zr,r(s) are generalizations of <f)(s) and Zr(s), respectively. Actually, one can see that 
4>i(s) = 4>(s), Zr,i(s) = Zr(s), «i,i(i) = —2 and 0:1,2 (i) = 4 (see Remark 3.10) and hence the 
equation (1.5) (resp. (1.6)) coincides with (1.1) (resp. (1.2)) when r = 1. 

We remark that, from the equation (1.5), remarkable cancellation of the singularities on (—00, — 1] 
occurs (though both functions Grp r ( s ) an d Zr r( s ) have singularities on the above interval (see Fig- 
ure 2 and 3), their product Dr r (s) does not). Visit also Remark 4.7 for some observations on these 
singularities. 

The organization of the paper is as follows. In Section 2, from the integral representation 
of the spectral zeta function, we first show the existence of Dr, r (s) and, then, study its basic 
analytic properties. In particular, we give a proof of the claim (i) in Theorem 1.1 (Theorem 2.3). 
Here, we make a special choice (2.4) of the log-branch which allows us reasonably to reach the 
functional equation of Dp* r (s) (see Remark 2.4). Next, using the Selberg trace formula, we derive 
the factorization (1.5). In Section 3, we determine the explicit expression of the gamma factor 
4> r (s) in terms of the Milnor gamma function T r {z) (Proposition 3.6). Furthermore, we give two 
other expressions of <?V(s); one by the Barnes multiple gamma functions T n (z) (Theorem 3.9), and 
another by the Vigneras multiple gamma functions G n (z) (Corollary 3.11). We notice that the 
latter two take advantage of observing analytic properties of 4> r (s). Using ladder relations of the 
multiple gamma and sine functions, we also have a functional equation of 4> r (s) (Theorem 3.13). In 
Section 4, via the expression (1.5), we study analytic properties of the Milnor-Selberg zeta function 
Zv,r{s) such as an analytic continuation and a functional equation (Theorem 4.1). Moreover, 
introducing a "poly-Selberg zeta function" Z^ n \s) by a certain Euler product (which is regarded 
as another generalization of the Selberg zeta function), we show that Zr, r (s) can be expressed as a 
product and quotient of Z^ 7l \s) (Theorem 4.11). This is nothing but the Euler product expression 
of Zr,r(s). 

In the course of the explicit determination of the gamma factor (p r (s), we will encounter the 
following series involving the Hurwitz zeta function; 

Rm(t,z) :=E c 9 (2 i +1 ;1 ^ +m+1 - 

2 j + m + 1 

This type of series has been studied in several places (see, e.g., [A, KS]). It will be shown that the 
exponential of such a series (with z = |) is expressed as a product of the Milnor gamma function 
(Proposition 3.4). We also note that a similar discussion developed in this paper yields explicit 
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expressions of the higher depth determinants of the Laplacian on the higher dimensional spheres 
([Y]). See also [WY] for number theoretic analogues of the present study. 

In this paper, as usual, C, M, Q are respectively denoted by the field of all complex, real and 
rational numbers. We also use the notations Z, N and No to denote the set of all rational, positive 
and non-negative integers, respectively. 

2 Higher depth determinants 

2.1 Existence and basic properties 

Let Ar be the Laplacian on the compact Riemann surface r\H of genus g > 2, = Ao < Ai < 

(T) 

A2 < • • • — > 00 the eigenvalues of Ap and Spec(Ar) := {Xj \j G No}. For T > 0, define 9 Ar (t) := 
X^A i >T e ~ Xjt - It is known that 9^(t) converges absolutely for Re (t) > and has the asymptotic 
formula 9^(t) ~ E^°=-i c " T) * n as * 4- with c -i ^ ( see > e -g-> [ R ])- Notice that c^} is not 
depend on T and hence we denote it by c_i. Define £X r (w,z) := + z)~ w where we let 

(Aj + z)~ w := exp(— w log {Xj + z)) for all j G No- In particular, we put CA r ( w ^ z ) := Ca^( w i z )- 
The series Q^{w, z) converges absolutely and uniformly in any compact set in {w G C | Re (w) > 
1} x {z G C I Re (z) > — T} and hence defines a holomorphic function in the region. At first, one 
easily sees the following 

(T) 

Lemma 2.1. The function Q'(w,z) admits a meromorphic continuation to the region Cx{z£ 
C j Re (z) > —T} with a simple pole at w = 1 and being regular otherwise. In other words, it can 
be written as 



where tp^ (w, z) is a holomorphic function in C x {z G C | Re (z) > —T}. 



Proof Let Re (w) > 1 and Re (z) > —T. Then, from the integral expression 



(2.1) 




for any N £ No, we have 






+ 



F(w) 



1 



TV 

n=-l Jo 
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whence 

N 



dt 

T 



(2-2) Ci>, z) = ^ / t»e-«(0$>(t) - £ <£°* 

. 1 y y cf } (-z)- i r tWe - z t 9 (T) (t) dt 

T(w) ^ ^ m!(w + n + m) iYio) A Ar K J t ' 

v ' n=-l m=0 v ' 

Since 6>2?(i) - ^ = _ 1 cl T) i n = O^ 4 - 1 ) as i 4- 0, the first integral in the righthand-side of (2.2) 
converges absolutely for Re (to) > — (N + 1) and hence, as a function of w, defines a holomorphic 
function in the region. The second term defines a meromorphic function on C having a simple pole 

(T) 

at w = 1 (notice that the points to = 0, — 1, —2, . . . are not poles of <^ r (it?, z) because of the gamma 
factor). Moreover, since Re (z) > —T, the last integral converges absolutely for all w G C, whence 
it defines an entire function. Therefore, letting N — > oo, we obtain a meromorphic continuation of 
CzjK z) to C x {z G C | Re (z) > — T}. □ 

We now study the higher depth determinants. Prom Lemma 2.1, the function 



Detf)(A r + .) := ex P (-|-Cf>, z) 



w=l—r 



is well-defined and is holomorphic in {z G C | Re (z) > — T 1 }. Notice that, since 

(T) 

is holomorphic, by the definition, Detr ( Ap + z) has no zeros in the region. In particular, we put 
Det r (A r + z) := De4 0) (A r + z). 

Proposition 2.2. The function Det r (Ar + z) can be continued analytically to 

(i) an entire function with zeros at z = —Xj when r = 1. 

(ii) a holomorphic function inC \ (— oo,0] with no zeros when r > 2. 
Proof. Let Re (z) > 0. Then, for any T > 0, we have 

Ca f (w,z) = E (Ay + ^ + CgW)- 

0<Aj<T 

From Lemma 2.1, the righthand-side is holomorphic at w = 1 — r. Hence, differentiating both sides 
at w = 1 - r, we have 

(2.3) Det r (A r + z) = JJ exp((Aj + zf" 1 log(A j + z)) • Detf )(A r + z). 

0<Xj<T 

When r = 1, since exp(log (Aj + z)) = Aj + z, the first factor in the righthand-side of (2.3) is a 
polynomial and hence is entire. This shows that Deti(Ar + z) can be extended analytically to 
{z G C | Re (z) > — T}. When r > 2, the first factor in this case is holomorphic in flo<A <r ^ \ 
(— oo, — Aj] = C \ (— oo, 0], whence Det r (Ar + z) can be extended to {z G C | Re (z) > — T} \ (-T, 0]. 
Therefore, letting T — > oo, one obtains the desired claim. □ 

Based on the above discussions, we next study the case z = — s (1 - s). Write Aj = r| + \ 
where rj G iM>o if < Aj < | and rj > otherwise. Moreover, let := ^ ± irj. We notice 
that Aj — s(l — s) = (s — ad~)(s — aj) and <x^ G [0,1] with < aj if < Aj < j and 
Re(a^) = | otherwise. As is the case of the previous discussion, we also start from the zeta 
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function Q A J(w, — s(l — s)) = ^A^T^i ~~ ~ s ))~ w > however, in this case, for small j, we replace 
the branch of the logarithm. Namely, we let (Xj — s(l — s))~ w := exp(— wlog^\Xj — s(l — s))) 
where 



(2.4) log^A,- - s(l - s)) := log |A,- - s(l 



+ . f (arg + (s - at) + arg_(s - 0j )) (0 < Xj < ±), 
\arg(( S -at)( s - a -)) (A, > ±) 

with the argument arg ± being respectively taken as — < arg + z <\ix and — |vr < arg_ z <\is 
for z G C. We notice that log^ = log if Xj > \. It is easy to see that lj(s) := log^(Aj — s(l — s)) 
is a single-valued holomorphic function in Wj where 




too 



) U [ 0j ,a j + ioo)) (0 < Xj < {), 



3 l / 1 1 

ioo) U [aj, - - ioo)J (Aj > |), 

and satisfies the trivial functional equation lj(l — s) = lj(s) for s G Wj. Let £/( T ) := {s G 
C|Re(— s(l — s)) > — T} (see Figures 4, 5 and 6 for T = 0, T = j and T > |, respectively). 
Note that has two connected components if < T < j and is connected otherwise. Clearly, 
limr^oo = C. Moreover, let := C\x >T W j- Notice that D U^ T \ From Lemma 2.1, 

the function (^(w, — s(l — s)) is meromorphic in C x (W^ fl U^) = C x and is in particular 
holomorphic at to = 1 — r for all r G N. Let Dp T ^(s) := Det£^(Ar — s(l — s)), which is holomorphic 
in and satisfies D^J(1 - a) = d{?J(*) for s G t/ (T) . Define the region O r := f]f =0 Wj. Notice 
that Or is equal to the righthand-side of (1.4) and remark that it is connected if and only if 
i G' Spec (Ap). Similarly to the previous discussion, an analytic continuation of Dr, r (s) := Dp°j,(s) 
beyond the region is given as follows. Let jr := min{j | Xj > \} and T > Xj F {> \). For 

s G U, we have 

(2.5) D r , r ( S ) = e r T r ) ( S )-D r i ; ) ( S ), 

where 

4?W-'= II exp((A J - S (l-^)) r " 1 log (j) (A,- S (l- S ))). 

0<Aj<T 

Note that, by the definition, e r (s) also satisfies the trivial functional equation e r r (1 — s) = e r r (s). 

When r = 1, since e r -[(s) is again a polynomial, it is an entire function. On the other hand, when 

r > 2, e rr (s) is holomorphic in Qr- From the equation (2.5), these show that Dr, r (s) can be 
extended to if r = 1 and n Or otherwise. Therefore, letting T — > oo and noting that 
limr^oo f/( T ) n Or = Or, one obtains the following 

Theorem 2.3. The function Dr, r (s) can be continued analytically to 

(i) an entire function with zeros at s = oc^ when r = 1. 

(ii) a holomorphic function in Or with no zeros when r > 2. 

Moreover, Dr, r (s) satisfies the functional equation Dr, r (l — s) = Dr r (s). □ 
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Remark 2.4. Define the region := U n {s G C | ± Re (s) > ±\}. Notice that U = U + U U~ . 
Since the function Dr r ( s ) is defined on each region U^, let us write the restriction of Dr r to 
as Dp . It is clear that 

(2.6) D*.(1-*) = D±.( 8 ) (se^). 

When r > 2 and | ^ Spec(Ar), the statements of Theorem 2.3 are equivalent to the following. 
Both functions D rr (s) can be continued analytically to the region Or and the equation Dp" r (s) = 
Dp r (s)(= Drj-(s)) (s G ^r) follows from the identity theorem. Hence, from (2.6), the functional 
equation of Dr, r (s) is immediate. 

Remark 2.5. Let r > 2. Then, the function Dr,r(s) is also extended beyond the region U even if 
we take the principal branch of the logarithm to define (Xj — s(l — s))~ w for all j G No- Actually, 
by the same manner as above, one can show that Dr^(s) admits an analytic continuation to the 
region C\([0,l]U(i+iR)), which never become connected for any V. Moreover, one finds that our 
choice of the log-branch is compatible with a functional equation and an Euler product expression 
of the Milnor-Selberg zeta function Zr, r (s) established in Section 4. 

2.2 Calculations of the spectral zeta function 

In this subsection, we study the function Ca. f (w, —s(1 — s)) by using the Selberg trace formula 
(2-7) f; f(rj) = £ lQ g ^7) r / (i ogiV ( 7 )) + ( 5 _ 1) /(r)rtanh(7rr)dr. 

,=0 7 6Hyp(r) ^(7) 5 -^(7) 5 J-°o 

Here, Hyp (T) is the set of all hyperbolic conjugacy classes in T, 5 7 G Prim (T) for 7 G Hyp (T) is the 
unique element satisfying 7 = 5^ for some k > 1 and / is a test function whose Fourier transform 
/( r ) := f(x)e~ trx dx satisfies the conditions / is holomorphic in the band {r G C | |Imr| < 
\ + f(—r) = f( r ) an d f(r) = 0(\r\~ 2 ~ 8 ) as \r\ — > 00 for some 6 > 0, 

Let Re (w) > r with r G N and s G U. Throughout the present paper, we always denote by 
t = t(s) := s — \. From (2.1) with T = 0, we have 

CA r {w + 1 - r, -s(l - ,)) = — -± r f +i- e -* (J2 e^A f . 

T(w + 1 - r) Jo ^ J £ 
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Applying the trace formula (2.7) with the test function fix) = )—^ e 4 ? (then f(r) = e~ r *) to 
the inner sum in the above integral and changing the order of the integrations, we have 

/oo 
(x 2 + t 2 )~ w+ xta,nh(nx)dx 
-oo 

= er, r (»,f) + 2(9-l)W r " 1 )^- 1 - f )j M+1 ( Wit ) i 



where 



with 



T(w + 1 - r) J £ 



1 >^ lOg N(Sj) t 2g (l°g ^(7)) 



r (£,i):=— f= V 

2 V^ ^tt,^ N(rf)% - N(-yyt 



7 eH yp (r) 
and 

/•oo 

(2.9) J m (w,t) : = / ip m (x;w,t)dx 

JO 

with 

ip m (x;w,t) := (x 2 + t 2 )~ w x m tanh(7rx). 

Here, (x 2 + t 2 ) -1 " := exp(— u; log (x 2 + t 2 )). We call J m (w,t) an m-th moment function. When 
s (z U, since Re(i 2 ) > j and there exists a constant er > such that log N("f) > er for all 
7 £ Hyp(r), the integral in Qr,r(w,t) converges absolutely for all w £ C and hence Qr, r {w,t) 
defines an entire function as a function of w. Moreover, we will show in the next section that 
J m (w, t) can be continued meromorphically to the whole w-plane C and is in particular holomorphic 
at w = 0. Therefore, differentiating both sides of (2.8) at w = 0, we see that the higher depth 
determinant Dr, r (s) can be expressed as 

(2.10) D r , r (s) = Gr,r(a)^,r(«), 
where Gr, r (s) '■= ^(s) 9-1 with 

Vt 1 ( 9 N 2( r 7 1 )* a Cr-i-<) 

(2.11) 4> r (s) := [[exp^- — J 2 e+i(w,t) 

£=o W 



w=0 



and 

(2.12) Zr, r (a) : = exp(--^-e r ,r(w, t) 

Our next task is to examine each factor i^ r (s) and Zr, r (s) more precisely. 



3 Gamma factors r (s) 

3.1 Derivative of the moment function J m (w,t) at w = 

To obtain an explicit expression of the gamma factor (j) r {s), let us study the function J m (w,t) 
defined by (2.9). Notice that, for a fixed s G C, the integral J m (w,t) converges absolutely for 
Re (w) > ^jr^ and hence, as a function of w, defines a holomorphic function in the region. In 
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what follows, from the definition (2.11), we may assume that m is odd. Moreover, for simplicity, 
we assume that ^ < s < 1, which implies that < t < \. Then, as a function of z, (p m (z',w,t) 
is holomorphic in {z £ C | Re (z) > 0} and is naturally continued meromorphically to the region 
C \ ((-oo,0] U [-it, it]) by extending log (z 2 + t 2 ). 

Suppose that Re (w) > Let R be a sufficiently large positive integer and e a small real 

number. Consider the counterclockwise contour integral with the integrand ip m (z;u),t) (on the 
z-plane) along the path 

C(t; R, e) := [+0, R] U C + (0; R) U [-R, -0] U L_ (t; e) U C(it; e) U L+(t; e), 

where C + (0;i?) is the semi-circle in the upper half plane with radius R centered at the origin, 
C(it; e) is the circle with radius e centered at it and L + (t; e) (resp. L_(t; e)) is the right (resp. left) 
side of the segment connecting and (1 — e)it. We take e so small that the circle C(it;e) contains 
no points of the form i(k + i), that is, the poles of tanh(-7rz), for all k E Nq (see Figure 7). 




1 4 




-R 


Figure 7 


: C(t;R,e). 


R 



Rc 



The residue theorem yields 

(3.1) / 
Jc 



R-l 



ip, 



>C(t;R,e) 

On the other hand, we have 



(z;w,t)dz = 2i m+1 e- mw Y,( k + -)' 2w+m (l -t 2 {k + -)~ : 



A;=0 



(3.2) / ip m (z;w,t)dz = / ip m (z;w,t)dz+[ / + / 

JC{t;R,e) JC+{0;R) \J[-R,-0] J [-+ 



'C+(0;R) \J[-R,-0] J[+0,R] i 

+ A m (w, t; e) + B m (w, t; e), 



p m (z;w,t)dz 



where 



A m (w,t;e) := / + / \<p m (z;w,t)dz, 

\JL + {t;e) JL-(t;e)J 

B m (w,t;e) := / ip m (z;w,t)dz. 

JC(it;e) 

Now, let us calculate each integral in the righthand-side of (3.2). At first, it is easy to see that the 
integral on C + (0;i?) converges to as R -> +oo because it is 0(-R m+1-3Re M). Next, since m is 
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odd and arg(z + t 2 ) = (resp. 2tt) if z G [+0, ii] (resp. 2 G [-R, —0]), we have 



|¥> m (z;ii;,i)dz = (1 + e 2 ™) / ip m (x;w,t)dx. 
[-R-0] J[+0,R]J JO 

r R 

= 2e mw cos(ttw) / tp m (x;w,t)dx. 
Jo 



This shows that 



lim / + / \ip m (z;w,t)dz = 2e mw cos(ttw) J m (w, t). 

R -^+°°\J[~R-0] J[+0,R]J 

Finally, we calculate A m (w,t;e) and B m (w,t;e). Notice that both functions are entire as functions 
of w and, by the Cauchy integral theorem, do not depend on the choice of e. Since arg(z 2 + t 2 ) = 
(resp. 2ir) if z G L + (t;e) (resp. z G L_(i;e)), we have 

/•(l-E)t 

t; e) = (1 - e" 2 ™") / y> m (i£; 

(l-e)t e2 







2i m+1 t- 2w e- wiw sm{7rw) (l - ^)~ w £ m tan(vr£K 



whence 



rt ti 

A m (w,t;e) = lim A m (w,t;s) = 2i m+1 r 2w e- mw sm{>Kw) / (l - \) w f" tan«)(i£. 

Moreover, by straightforward calculations, we see that B m (w,t;e) = 0(e 1 - Re ^)) as e — > 0. This 
shows that B m (w,t;e) = lim e _>o B m (w, t; e) = for Re (w) < 1, whence, by the identity theorem, 
B m (w, t;e) =0 for all w G C. Therefore, letting — > +oo in the formulas (3.1) and (3.2), we have 

jm+l . . 

(3.3) Jm(w,t) = -(J m l(w,t) + J m ,2(w,t) , 

COS 7TW V / 



where 

(3.4) J mA (w, t) := -t~ 2w S m(Trw) [ (l - J) ~™£ m tan^)^, 

jo 4 

(3.5) J m , 2 (w,t) := J> + I)^ +m (l - t 2 (k + i)- 2 )^. 

fc=0 

The function J m i(w,t) is clearly entire as a function of w. Let us calculate the derivative of 
Jm,i(w,t) at w = 0. To do that, we employ the following basic multiple trigonometric functions 

9 n 

([KOW1], see also [KK]). Put Py(u) := (1 - u) and P n {u) := (1 - u)exp(u + \ + • • • + \) for 
n > 2. Then, the basic multiple sine S n (z) and cosine function C n (z) are respectively defined by 

mSZ 

<S n (z) : 



2 ^ n = 2 ^ iK 1 - ^) = 2sin ^) ( n = *)> 

mSZ m=l 

M^d II ) ( - ir T~' <" > 2). 

\n — 1/ - LJ - V v m v m / 



mSZ 
2 \ ("2") 



W:= J] 2 =5 n (2,) 2l ^5 n (^)- 1 . 



m:odd 
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Notice that, when n > 2, S n (0) = C n (0) = 1. Moreover, for z / 0, we have the following integral 
expressions. 

(3.6) S n {z) = exp(^ < n ~ X cot«)de) , C n (z) = exp(- ^ 7rf n_1 tan«)cZ£ 
Proposition 3.1. It holds that 

d 

(3.7) -z-J m ,i(w,t) =logC m+ i(i). 
Proof. From (3.4), we have 

J m ,i(w,t) = -(1 + 0(w))(ttw + 0(w 3 )) [ (l + OH)e m tan«)de 

Jo 

7rC m tan«)^ )w + 0(u> 2 ). 



Hence, the claim follows from the formula (3.6). □ 

We next study the function J mj 2(w,t). We first show the following 
Lemma 3.2. We have 

(3.8) J m>2 (w, t) = f2( W + J ~ r ) t 2j ({2w + 2j -mi). 

i=o V 3 ' 

This gives a meromorphic continuation to the whole plane C as a function of w with possible simple 
poles at w = — j for j G No- In particular, it is holomorphic at w = 0. 

Proof. Notice that, since < t < \ (recall that we assume that ^ < s < 1), it holds that 
\t 2 (k + l)~ 2 \ < 1 for all k > 0. Therefore, from (3.5), using the binomial theorem, we have 

= E(* + 5)" 2 " + ™ £ (" + ■ " 1 W + \r'" 

k=0 j=0 \ J / 

-E(" + r i )«"<(*»+!«— .5)- 

Hence one obtains the expression (3.8). Moreover, since ((2w + 2j — m, ^) is uniformly bounded 
with respect to j, this gives a meromorphic continuation to C. Now the rest of assertions is clear 
because the Hurwitz zeta function C( w , z ) has a simple pole at w = 1. Notice that the origin is not 
a pole of Jm,2{w,t) since, when j = > 1, ( w+ j _1 ) = 0(w) as io — > 0. □ 

Before calculating the derivative of J m ,2(w,t) at w = 0, let us recall the multiple gamma 
functions, which is defined via the Barnes multiple zeta function ([B2]) 

(n(w,z) := — (Refit;) > n). 

mi,...,m„>0 

This clearly gives a generalization of the Hurwitz zeta function; £1(10, z) = £(w,z). It is known 
that Cn( w , z ) can be continued meromorphically to the whole plane C with possible simple poles at 
w = 1, 2, . . . , n. The multiple gamma function T n<r (z) of depth r is defined by 





( v 



w=l—r 
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In particular, we put T n (z) := T nt i(z) and T r (z) := Fi ;r (z). These are respectably called the Barnes 
multiple gamma function ([B2]) and the Milnor gamma function of depth r ([M], see also [K0W2]). 
From the Lerch formula -§^C,( w i z )\ w -o = 1°& we nave ^l,l( z ) = T±(z) = Ti(z) = whence 
these in fact give generalizations of the classical gamma function. We remark that T n (z)~ l is an 
entire function with zeros at z = — k of order (^"Y 1 ) f° r k € No. 
From the expression (3.8), it can be written as 



m — 1 
2 



1 f W I l\ 1 

■WtM) =<Z(2w-m,-) + ( 3 ~ )t 2 K{lw + 2j-m,-) 

3=1 

+ (" % " > m+1 f 2l ° + »• 5) + £ (" + j " + 2; - m, \) 

=: T!(w, t) + T 2 (w, t) + T 3 (w, t) + T 4 {w, t). 



Using the expansions { w+ j X ) = j(w + H(J — l)w 2 + 0(w 3 )) as u; — > for j > 1 where H{m) : = 
YlT=i \ an d CC 1 ", = - ^(z) + C(w - 1) as to ->• 1 where ip(z) := ^logr(z) = jr( z ) is the 
digamma function and employing the formula Q(l — m, z) = — Bm ^f > for m € N where B m (z) is the 
Bernoulli polynomial defined by ^§~zf = Ylm=o ^m{ z )^-\i we have 

(3.9) T t (w,t) =C(-m i) + A C ( 2u; _ m ,I) «; + 0(w; 2 ) 

2 aw 2 7 - 



io=0 



i?m+l( ^ ) + (21ogr m+1 (i))u; + 0(u; 2 ), 



m + 1 V 2 



m — l 
2 



(3.10) r 2 («;, i) = -.{w + 0(w 2 ))t 2 ^ (C(2j - m, i) + 0(«/ 



i=i J 



j^j{m+l-2j) ) 1 h 



im+X 



771 — ^ 



(3.11) r s («;, t) = (w + H(^)w 2 + 0(w 3 )) f— - ^(5) + 0(^ 2 ) 



2 ; v VV2io v 2 

1 _ t m+i + / 2 n~,m-U 



r(^( ! V)-^^)) tm+1 ) u; + 0(u;2) ' 



m + 1 \m + 

00 .. -. 
(3.12) T 4 («;,t)= J] -{w + 0(w 2 ))t 2 ^C{2j-m,-)+0(w 

■> 2 

= (2R m {t))w + 0(w 2 ) 
as t/j — > 0. Here, R m (t) is defined by 



CO 



R 7*\ ._ V £( 2j * + 1 ' 2) f 2 7 +m+l 



This yields the following 
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Proposition 3.3. It holds that 



(3.13) ^Lj m2 (w,t) 



w=0 



m+l / \ -, 

- logC m+ i(t) + z^{ H C^) ~ 2H(m))t m+1 . 



m + 1 V v 2 

Proo/. This follows from the identity ^ J m ,2(w, t) \ w=0 = Ylj=i m T i( w ' t )\w=a together with (3.9), 
(3.10), (3.11), (3.12) and the following proposition, which will be proved in Subsection 3.4. □ 



Proposition 3.4. Let m> 1 be an odd integer. Then, we have 

m+l , x 11 

(3.14) R m (t) = £ ™ t m+l ~ k log T k (t + i) - - log C 

fc=i V" / 



(t) 



*ri j(m + 1 - 2j) 

Prom the equation (3.3), noting that cos(-7rw)~ 1 = 1 + 0(w 2 ) as w — > 0, we have 



^J m (u;,t) 



= i m+1 (l^J m ,l(w,t) + Wm, 2 (w,t) 

w =o \ OW w=0 OW 



w=0 



Therefore, substituting (3.7) and (3.13) into this equation, one eventually obtains the derivative of 
the moment function J m (w,t) at w = 0. 



Proposition 3.5. It holds that 



(3.15) 



•m+l 



, m — 1 , 



u;=o m + l\ 2 

m+l 



2H{m))t 



/.m+l 



/ \ 1 

k=i ^ ' 



□ 



3.2 Explicit expressions of r (s) 

We obtain the following expression of (p r (s) by the Milnor gamma functions. 
Proposition 3.6. We have 

2r 



(3.16) 



( S ) = e -^B^.^ 2r jj r fc ( fl )(^) f (-D^- 1 ^-!) 



r \ 2k f_i\k+r-lfr,„_i\2r-k 



k=r 



To prove this, we need the following lemmas about sums of the binomial coefficients. 
Lemma 3.7. Let r £ N. Then, the following equality holds; 

(2r)H 



(3.17) a: = E(% ^TJlW) " 2#(2* + 1)) 



r 2 (2r- 1)!!' 
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Proof. Since H(£) - 2H(2£ + 1) = -2 Yf k =o 2FFI ( we understand that H(£) = if I < 0), changing 
the order of the summations, we have 



(3.18) 



1 yVr-l\ (-!)*_ 2^/r-l\ 
'^2^ + 1^1 £ /^ + 1 r^l / 

i—n t—u V / i— n V / 



jfe 7 2fc + l 



'r-l\(-l) £ _ 2 ^ - 1\ (-l) fc 
Here, we have used the formula Yl^Zk C*^ 1 ) Tm 



,r 1 C ^ - ^ — - p ( r fc 1 ) ' which is easily obtained by induction 
on fc. Therefore, since the sum on the rightmost-hand side of (3.18) is equal to the beta integral 



\i - jy-ite = W- r) = r( ^ )r(r) = (2r)!! 

/ ' 2^ l 2' rj 2r(r + i) 2r(2r-l)!!' 



one obtains the desired formula. 



□ 



Lemma 3.8. (i) Let r G N and 1 < k < 2r. Then, the following equality holds; 

r-l 



(3.19) 



r \2k 
k — r I r 



— / j\fe+r— 1q2t— k 



In particular, D r (k) = 0ifl<k<r — 1. 

(ii) Let 1 < p < 2r. Then, the following equality holds; 



2r 



(3.20) 



k=p 



k — p 







(p : odd), 



Proof. From the binomial expansion 

4(_l)* ^-i / _j> 



(3.21) 



D r (k) 



(k-l)\dt k 

4(-l) fc 
(fc = l)! 



t=i 

jk-1 ik-2 



f=l 



This shows that D r {k) = 0ifl<fc<r — 1. Now, let r < < 2r and j > r — 1. Then, by the 
Leibniz rule, we have 



(l-t 2 ) 



2\r— 1 



(l-tr^i+t) 



r-l 



t=l 



E 

m=0 



mj \dt 
r — 1 



j - r + 1 



-(1-tf 



r-l 



t=l Vefc? 



r-l 



t=l 



.^\r-l22r-2-j 



Hence, from (3.21), using this formula with j = k — 1 and k — 2, one obtains (3.19). Moreover, it is 
clear that the equation (3.19) is also valid for 1 < k < r — 1 because ) = for such a k. Hence 
the claim (i) follows. We next show the claim (ii). Consider the generating function Ylp=i D r ,pX p . 
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Changing the order of the summations and using the formula (3.19), we see that this is equal to 

2r 2r lX 2r 



EEL- J D ^ xP = E a-(*mi + ^) fe ~ 1 
p =i fc= P ^ ^ fe=i 

= 4x 2 (l-x 2 ) r - 1 



'SO 15 " 



Hence the claim follows. □ 

We now give the proof of Proposition 3.6. 
Proof of Proposition 3.6. From the equation (3.15), changing the order of the products, we have 



(3.22) Ms) = <p* r nn r >(t+ y-^KKisv 



e=o k=i 

D C r t 2r TT -p / \D r (k)t 2r ~ k 



k=l 

Therefore, one immediately obtains the formula from (3.17) and (3.19). □ 

To clarify the analytic properties of 4> r (s), let us rewrite the expression (3.16) in terms of the 
Barnes multiple gamma functions. The following expression is obtained in [KOW2]; 

r 

(3.23) T r (z) = n r iW Vj ' W . 
where, for r > 1 and j > 1, c r> j(z) is the polynomial in z defined by 

(3.24) crj(z) := ( J ~i l ) " 1 ~ VT 1 - 

1=0 ^ ' 

For example, c rjT .(z) = (r — 1)!, cv^-^z) = \{2z — r)(r — 1)!, . . . and c r> i(z) = (z — l) r_1 . It is easy 
to see that the polynomial c r j{z) satisfies the recursion formula 

c r,j{ z ) = {z~ l)Cr-lj(^) + (j - l)Cr-l,j-l(z - 1). 

From this, noting that c±j(z) = (1 — l)- 7-1 = for j > 2, we have c r j(z) = for 1 < r < j — 1. 
Notice that c r j (z) is also given by the generating function 



(3.25) (T + zf' 1 =Y J C r ,j{z) 



r + i-i 
i-i 



For let us denote [xj by the largest integer not exceeding x. Then, we have the following 

expression of 4> r (s) in terms of the Barnes multiple gamma functions. 
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Theorem 3.9. We have 



2r 

(3.26) Ms) = e~^'-U (s ~^ )2r flTjis)^^ 3 -^, 

3=1 

where a r j(t) is the even polynomial defined by 
(3-27) a rd (t): = 4 £ " J) (-l)^W^ 

This gives 

(i) a meromorphic continuation of <f) r (s) to the whole plane C with poles at s = —k of order 
2(2k + 1) for fceN when r = 1. 

(ii) an analytic continuation to the region C \ ((— oo, — 1] U [0, —zoo)) when r > 2. 
Proof Prom the equations (3.22) and (3.23), it suffices to show that 

2r 

(3.28) a rJ (t) = Y J CkA t + ^)D r (k)t 2r ~ k . 

k=j 

Actually, from the equation (3.24), the righthand-side of (3.28) is rewritten as 

2r /j-l 



2r-2£ 



- e (' 7 *) (-1)' (e e C; ') (-5 - if-^D^-v'-A . 

Z=0 V 7 \k=j m=0 V 7 / 

Putting — m = p in the inner sums, one sees that this is equal to 

e e (tiMlEfVW^-H^ 



.k — p J \ ^— ' V I 

p=l \fc=max{pj} / \i=0 



e(_e 



P— 1 \fc=max{p,j'} 

Moreover, the summand for p = 1, 2, . . . , j — 1 vanishes because c p j(z) = 0, whence, consequently, 
this can be written as J2p=j D r ,pCp,j{\)t 2r ~ p ■ Therefore, from (3.20), one obtains the equation 
(3.28). The rest of the assertion follows from the equations 0:1,1 (i) = —2 and 0:1,2 (i) = 4 and the 
fact that the point z = —k is a zero of Tj{z)~ l of order ( k ~^^ 1 ). This ends the proof. □ 

Remark 3.10. From the equations c r , r „i(z) = \{2z — r)(r — 1)! and c r , r (z) = (r — 1)!, one sees 
that both a r flr-i(t) and av,2r(£) are integers respectively given by 

o r , 2r _i(t) = 4(-l) T - 1 c 2r>2r _i(i) = (-l) r (4r - 2)(2r - 1)!, 
a r , 2r {t) = 4{-l) r c 2r , 2r (\) = 4(-l)'- 1 (2r - 1)!. 
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We next give an expression of (f) r (s) via the Vigneras multiple gamma function G n (z), which are 
characterized by a generalization of the Bohr-Mollerup theorem ([Vi]). Notice that G\{z) = T(z) 
and G 2 (.z) = G{z) where G(z) is the Barnes G-function studied in [Bl]. From [SC, p. 87 (27)], we 
know that G n (z) is essentially equal to F n (z); 

(3.29) G n (z) = e (- 1 )" E"r H nj ( 2 )C(- J ) . r^-i)- 1 . 

Here, b n j(z) is the polynomial of degree n — 1 — j defined by the generating function ( 3 '^™7 ) = 
Y^k=Q^n,k{z)(j + z) k . To be more precise, let s(n,m) be the Stirling number of the first kind 
defined by (z) n = ^(-lj^sfn, m)z m where {z) n := = z(z + 1) • • • (z + n - 1) is the 

Pochhammer symbol. Then, it is given by b n ^{z) = ^-j^m ^I=t Cfc) s ( n ' m + ^-)z m ~ k ■ The 
following expression is immediately obtained from (3.26) together with the identity (3.29). 



Corollary 3.11. We have 



2r 

(3.30) Ws) = e -^( s 4)^^(HK'H) nGj(s)( . riarj(s4)i 

Here, I3 r ^{t) is the polynomial defined by 



2r 1 



3=1+1 



□ 



Example 3.12. Let t = s - \. Then, from (3.26) and (3.30), we have 



Ms) = e- 2t r!( s )- 2 r 2 ( s ) 4 

= e -2^-4tC'(0)+4C'(-l) Gi(s) -2 G2(sr 4 ) 

Ms) = e -l i4 r 1 ( s )^* 2+ 5r 2 ( s ) 4 * 2 - 13 r 3 ( s ) 36 r 4 ( s )- 24 

= e -|t4-8^C'(-l)+12*C'(-2)-4CX-3) Gl(s) -2^+| G2(s) -4tHl3 G3(s) 36 G4(s) 24 ) 

3 ( s ) = e-I^r!^)- 2 ^ 2 -^*) 4 * 4 - 2 ^ 



= e -i^-16t 3 C'(-2)+32t 2 C'(-3)-20tC'(-4)+4C'(-5) 

xG 1 ( s )- 2t4+i2 -lG 2 ( s )- 4 * 4+26 * 2 -^G3( S ) 72 ' 2 - 330 G4( S ) 48 ' 2 - 1020 G 5 ( S )-^ 
Notice that the expression of </>(s) = </>i(s) is obtained in [Vo]. 

3.3 Functional equations of (p r (s) 

In this subsection, we establish a functional equation of cp r (s). To do that, we first recall the 
normalized multiple sine function S n (z) studied in [KK]; 

(3.31) S n (z) := Tnizy'T^n - z)^" . 

Notice that, from the reflection formula, we have S\(z) = 2sin(-7rz) = S\(z). We remark that it is 
shown in [KK, Theorem 2.14] that S n (z) can be expressed as a product of Sj(z) for j = 1, 2, . . . , n 
and vice versa. 
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Theorem 3.13. We have 

2r 

(3.32) r (l - s) = (J] S k (s) a ^-^ Ms ) 

k=l 

for all 

(i) s£C when r = 1. 

(ii) s£C\ ((-oo, -1] U [2, +oo) U [0, -ioo) U [1, +ioo)) when r > 2. 
To obtain this, we need the following lemmas. 

Lemma 3.14. We have 

(3.33) r n (l -z) = fliSjiz^iz))^^. 
Proof. Let E n (z) := 5 n (z)r n (^). Note that, from the ladder relations 

(3.34) r n (z + i) = r n (z)r n _ 1 (z)- 1 , 

(3.35) S n (z + l) = S n {z)S n . l (zY\ 

we have -E^z + 1) = E n (z)E n -\(z) . Here, we put Tq{z) := z _1 and £0(2) := — 1. By the 
definition (3.31), we have S n (z)( -1 ^" = T n (n — z). Hence, using the relation (3.34) repeatedly, we 
have 

n-2 

E n (z)^ n = T n (n -z)= r n (l - z) XI r n _!(n - 1 - (* + m))" 1 

m=0 

n-2 

= T n {i-z)\[E n _ 1 {z + m)^ n . 

m=0 

Therefore, using the equation E n -i(z + m) = llj=n-i-m Ej(z) Kn-i-jJ^ which is obtained 

from the ladder relation of E n (z) and the formula X^m=a (T) = (b-a)> one sees th & t 

.!)» \ 



Y n {l-z) = E n (z)^ 1)n \ l[E n ^(z + 

\m=0 

/n-2 n-1 

=^(*) { - i)n (n n 



m=0 / 
n-2 n-1 \ -1 



V I ^—.m— n — 1— j \n—l—j/ 



\m=0 j'=n— 1— m 
n-1 

E n {z)^ l)n Y[ Ej(zV 

3=1 



= l[E j (zr - 

This shows the claim. □ 
Lemma 3.15. For 1 < k < 2r, we have 

(3-36) ar, k (t) = (-l) k J2 (I _ ijorM 

j=k V 7 
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Proof. Write the righthand-side of (3.36) as a r ^{t). It suffices to show that Y^k=i a r,k(f)^^-\ *) 
X]fcLi S r ,fc(0( T j^]" 1 )- 111 fact, from (3.27), using the formula (3.25), we have 

2r /T i „• 1 \ i /T 1 _i_ I\2 „ r _i 



2 ' 



On the other hand, using the identity Efc=i( _1 ) fe (lli) Tfc-i ^ = _ (~^-i~ 2 )> we have 



fc=l v " ' j=l 



1 / f-T - ±") 2 \ r- 1 



2 ; V *2 
4t—(T + i)fl ( r+ 5)V- ] 



2 7 V i 2 

Hence the claim follows. □ 

We now give a proof of Theorem 3.13. 

Proof of Theorem 3.13. We first notice that both functions t 2r and a r j(t) are invariant under the 
transform s t— > 1 — s because they are even polynomials in t. Hence, replacing s with 1 — s in (3.26) 
and using the formulas (3.33) and (3.36), we have 

r (i_ s ) = e ^r^ f [J m(5 fe ( s )r fc ( s )) ( H-J 

i=i Vfc=i / 



2r 

e -^^ f2 ''n(5 fc ( s )r fc( ,)) ( - 1)fc ^(-)^w 



k=l 



2r 

(n^( s ) ar ' fc(i) )^ (s) - 



k=i 

This shows the equation (3.32). □ 
3.4 Proof of Proposition 3.4 

Let m € Nq. The aim of this subsection is to give a proof of Proposition 3.4. Let 



2 7 + m + 1 



Note that R m {t) = R m (t, A). We start from the identity ([SC, p. 159 (4)]) 

m, z) = f; ^±kg) t2j+i = i ( log r(z _ t) _ log r(z + ^ + < ]^|). 
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Letting z = \, we have 

(3.37) i^(f) = ^V^A)(^H 



-* m (t, i) - ^logC m+1 (t) + -L_^)t m+1 , 



where 



<mm) := /Vv^+^e 

Jo 



Notice that, in the last equality, we have used the formula — £) = ^(^ + — 7rtan(7r£) and 
(3.6). Hence, it is enough to evaluate the integral <$> m (t,z). Define the polynomials n A_k(z) for 
1 < k < n — 1 and n B m (z) for m 6 No by the generating function 



+ Jn-z)t n ~ 1 00 /m 

- m B m (z)—.. 



V ; k=l m=0 

These are called the Barnes multiple Bernoulli polynomials ([B2]). Notice that the degree of n B m (z) 
is m + n — 1. Since \B m {z) = B m (z), n B m (z) gives a generalization of the Bernoulli polynomial. 
In fact, using the polynomial n B m (z), one can evaluate the special values of the Barnes multiple 
zeta function £ n (w,z) at non-positive integer points; 

(3.38) Ul-m,z) = - nBm ^ (meN). 

m 

To obtain an explicit expression of 3> m (i, z), we first show the following 
Lemma 3.16. It holds that 

(3.39) (m + 1) /VlogrV(£ + 2K 

JO 

= E X ) ( r + I " ^ "V* 1 "* l °sr n , r+k (t + z) + P n , r (t, z; m), 

where P n , r (t, z; m) is the polynomial in t of degree n + r + m defined by 

(3.40, PnA t, „ m) : _ £ t}U^±l2} (™ + + 'r^l, + ,) 

smj-i An + r\ 1 / „ , . H(r,m + r) „ . , 

+ ^ 1 (io g r njm+r+ i(* - 1 

with H(m, n) := Ylk=m k = ~~ H(m — 1) /or n > m. 

Proof. Integration by parts yields 

/ r(e+a)-^ = E / i\ T;, ? ( Iff — + (-i) m+1 -' 

Jo f— ' (m + 1 - £;)! (-w + 



^ (m + 1 - fc)! (-w + 1)* ' (~w + l) m+ i 
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Hence, for Re (w) > m + n + 1, changing the order of the integral and summation, we have 



mi,...,m n >0 
m+1 



E 



(_l)*-i m ! t^+^Cnfa-M + z) ^ m+1 ,C„(w-m-l,z) 



+ (-l) m+i m!- 



(-W + l)m+l 



^(m + l-fc)! (-w + l) fe 

This gives a meromorphic continuation of the lefthand-side to the whole plane C. Now the desired 
formula (3.39) immediately follows from the above equation by differentiating both sides at w = 1—r 
together with the following equation obtained from (3.38); 

ow\ {— s + ljfc / w=i-r (r)k V r + k 

Corollary 3.17. It holds that 

m+1 , ^ 

(3.41) $ m (t, z) = Y J ("l) fc+1 ( i.^ i ) lQ g r *(* + z ) + P -(*> z )> 

fc=i V * / 



□ 



where P m (t, z) is the polynomial in t of degree m+1 defined by 

m+1 , ^ {-^B E}±1 _ l {z), 



(3.42) P rn (t, z) : = —^—H (m)t m+1 + V 

m+1 



i< + i og r m+1 (z). 



^ /(m + 1-/) 
Proof. From the formula (3.39), we have 

<S> m (t,z)=t m logT 1 (t + z)-m /V^logrV^ + zK 

J o 

m+1 / \ 

= ^ (-l) fc+1 ™ i"* 1 "* logr fc (i + *) - P M (t, z; m - 1). 
fc=i V * / 

Here, we have used the identity i\i(z) = Ti(z) = Hence it suffices to show that 



(3.43) 



-P ljl (t,z;m-l) = P m (t,z). 



From the definition (3.40), using the identity B m (t + z) = YT=o i^)Bi{z)t m - 1 and changing the 
order of the summations, we have 



-P lt i(t,z;m - 1) 



(-l^^Vm + l 



m 



+ 



1=0 

(-l) m H{m) 
m + 1 



di{m)B m+ i_i(z)t l 



B m+l (z) + (-ir +1 io g r m+1 (z 



where di(m) := Ylk=o (jD( — l) k H(m — k). This implies that, to obtain the equation (3.43), it is 
enough to show that 



di(m) = < 



H(m) 



(1 = 0), 



^(m-QKJ-Dl 



ffi! 



[(-l) m+1 ^(m) 



(l=m + 1). 



24 



Nobushige KUROKAWA, Masato WAKAYAMA and Yoshinori YAMASAKI 



Actually, the case I = is clear. Let 1 < I < m + 1 and Fi(s) := Ylm=odi( m ) sTn a generating 
function of di(m). Then, from the identity Ylm=i H(m)s m = — (1 — log (1 — s), we have 

1 /A 00 

fc=0 ^ 7 m=l 

Hence, by the Leibniz rule, we have 



1 rt m 1 J^L /m\ rt m ~ k rl k 

fc=0 v 7 



s=0 



A; \m — k 

fc=max{m+l — 1,1} 



When 1 < I < m, this is equal to 

y ( l - *) = f 1 x ™-i (1 _ x -iy-i dx = -w-iy- 

f^ Q m-k\ k J J ml 



On the other hand, when I = m + 1, using the formula ip(m) = —7 + H(m — 1) for m £ N and the 
equation (see, e.g., [SC, p. 15 (13)]) 

1 _ 1 

_ de = _ 7 _^) (Re(z)>0), 



£ 

where 7 = 0.57721 ... is the Euler constant, we see that the rightmost-hand side of (3.44) equals 

( _i r vizl}!(Y) =(-ir C {1 ~ x)m ~ 1 dx = (-ir +i H( m ). 

k=1 k \kj J Q x 

This completes the proof. □ 

We now give the proof of Proposition 3.4. 

Proof of Proposition 3.4- Notice that, since m is odd, B m+ i_i(^) = if Z is odd. Hence one obtains 
the formula (3.14) from (3.37) together with the formulas (3.41) and (3.42) with z = ^- d 

4 Milnor-Selberg zeta functions Zp r (s) 

In this final section, we study the Milnor-Selberg zeta function Zy, t {s)- Recall that Zt^(s) is 
defined by (2.12) and is holomorphic in the region U (see Figure 4). 

4.1 Analytic properties of Zr, r (s) 

The following theorem is easily obtained from the earlier discussion. 

Theorem 4.1. The function Zp )7 .(s) can be continued analytically to 

(i) the whole plane C when r = 1. 

(ii) Or \ (—00, —1] when r > 2. 
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Moreover, it satisfies the functional equation 

2r -l 

(4.1) Z r>r (l " *) = (II Sj(s)- a ^ s -^Y Z r , r (s) 

3=1 

for all 

(i) s € C when r = 1. 

(ii) s£fi r \ ((-oo, -1] U [2, +oo)) w/ien r > 2. 
Proof. From (2.10), we have 

(4.2) Z r , r (s) = ^(s)^^" 1 )D r , r ( S ). 

This gives the desired analytic continuation from Theorem 2.3 and Theorem 3.9. The functional 
equation (4.1) immediately follows from Theorem 2.3 and Theorem 3.13. □ 

To study a "complete Milnor-Selberg zeta function", we need the following lemma. 

Lemma 4.2. For < I < 2r - 1, let 

2r-l /n 

2r - j 



&,,,(*) := (-1)' £ ( , •')orj(*). 



Then, for 1 < m < 2r, we have 

2r-l 



(4-3) E ( 2r ijM*) = (-l) m «r.m(*). 



Proof. Using the equation J2i= a (~ l)'(a) (/) = ^a.&( — l) a where <5 a) 6 is the Kronecker delta, we have 

2r— 1 / ; \ m / 2? — j , \ /o -\ \ 

E U 1 ra )M«) -EE (-D'U i J ( V) M*> = <- 1 '"' a -< ( »' 



/=2r— m j'=l \(=2r— m 



□ 



The following gives a generalization of the functional equation (1.3). 
Corollary 4.3. Define the complete Milnor-Selberg zeta function by 



2r-l 1 

9-1 



(4.4) 3r, r (a) := ( II ^ + O^"^) 2r,r(*)- 



1=0 

Then, we have 

(2r)H ( n iw„ l\2r 

(4.5) H r , r (s) =e*V» = W {8 K 2 > ■ D r , r (s). 

/n particular, Hr r ( s ) * s 

(i) an entire function when r = 1. 

(ii) a holomorphic function in Qr w/ien r > 2. 
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Moreover, it satisfies the functional equation 

Hr,r(l ~ s ) = s r,r(s) 

for all 

(i) s £ C when r = 1. 

(ii) s G Op w/ien r > 2. 

Proof. It is sufficient to prove the equation (4.5). To do that, from (2.10) with (3.26), it is enough 
to show that T 2r (s + l) & rM = J]^ =1 T m (s) ar ' m ^. Actually, from the ladder relation (3.34), 

one can show that F2 r (s + I) = Ylj=o ^2r-j(s) J ■ This yields 

2r-l 2r-l i 



1=0 1=0 j=0 



2r 



nr. 

m=l 
2r 

n r. 



m=l 
2r 



m=l 



In the last equality, we have used the equation (4.3). This completes the proof. □ 
Example 4.4. Let t = s — i. Then, we have 

E r ,i(*) = (r 2 (s) 2 r 2 ( S + lf) a - x Zr tl {s) = HrW, 

Sr , 2 ( s ) = (r 4 ( s )-i +2 * 2 r 4 ( s + i)-f - 2 * 2 r 4 ( s + 2)-" " 2 ' 2 r 4 ( s + sr^y-'z^s), 
e t>3 (s) = (r 6 ( s )l-* 2+2 * 4 r 6 ( s + i) 2 -¥-™ 2 -^r 6 ( s + 2 ) 2 t 1 + 22 * 2 + 4 ' 4 

x r 6 ( s + 3)^ +22i2+4i4 r 6 (, + 4)T- 2 " 2 -6* 4 r 6 ( s + 5)i-* 2+2 * 4 ) 9 " 1 z ri3 (,). 

Remark 4.5. Using the ladder relations (3.34) and (3.35), one can prove that 



a r-, j (t) 



j(s) rj — 1 - r 2r— 1 



This reads the definition (4.4) of the complete Milnor-Selberg zeta function. 

Remark 4.6. The equation (4.5) implies that we can get rid of the singularities of Zr, r (s) in 
(— oo, —1] by multiplying suitable gamma factors. Moreover, it also says that the Milnor-Selberg 
zeta function Zr, r (s) has no "non-trivial zeros" because Dr, r (s) does have no zeros. 

Remark 4.7. Let f(s) be a function on C and m(s) a polynomial. Suppose that it can be written 
as f(s) = (s — a) m ( s ^g(s) around s = a where g(s) is a holomorphic function at s = a with 
g(a) 0. In this case, let us say that /(s) has a "multiplicity polynomial m(s) at s = a" and write 
m(s) = m(s; f,a). It is clear that this is a generalization of the multiplicity (or order) of zeros or 
poles of meromorphic functions (they are the case m(s) £ Z). For example, if f(s) has a zero (resp. 
a pole) of order n at s = a, then the multiplicity polynomial of f(s) g ^ for a polynomial q(s) at 
s = a is given by m(s;f q ,a) = nq(s) (resp. —nq(s)). 
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From the expression (4.2), for k G N, one can see that 
(4.6) m{s; Z r , r , —k) = 2{g - l)(2Jfe + l)(s - k) r ~ l {s -k- l) r_1 . 

In particular, when r = 1, this coincides with the multiplicity 2(g — l)(2k + 1) of the trivial zero 
s = —k of the Selberg zeta function Z-p(s). The equation (4.6) is obtained as follows; from the 
equation (4.2) and the expression (3.26) together with the fact that Fj(s) -1 has a zero at s = —k 
of order ( fe ^7^), ^ can ^ e written as 

where zr, r {s) is some holomorphic function at s = —k with zr, r (— k) ^ 0. This shows that 

m(s; Zr>, -fc) = (g - 1) ( • _-, )«r,j(*) 

i=i V J ' 
2r 



4(5 - 1) E 



!=1 



r-1 



r I f_1 ^-l+2r-2€ 



2( 5 -i )( 2fc-i)^- 2 x;( r /) (-(^+^) 2 ^ 



£=0 

2 , t , l^- 1 



= 2(5- l)(2fc - l)^ 2 - (k + - 
Note that, in the third equality, we have employed the formula (3.25). 

4.2 Poly-Selberg zeta functions 

In order to study an Euler product expression of the Milnor-Selberg zeta function r (s), we now 
introduce a certain generalization of the Selberg zeta function. 

For m G N, define the function (s) by the following Euler product. 

oo 

(4.7) Z^\s):= UM N ( P r S - n f° SNiP)) ~ m+1 > 

PePrim(r) n=0 

where H m {z) := exp(— Li m (z)) with Li m (z) := ^ZfcLi jfer being the polylogarithm of degree m. 
This infinite product converges absolutely for Re(s) > 1. We call Z^ n \s) a poly-Selberg zeta 
function of degree m. Notice that, since Li\(z) = — log(l — z) and hence H\(z) = 1 — z, we have 
Z^\s) = Zr(s). To give an analytic continuation of Z < f l \s), we first show the following 

Lemma 4.8. It holds that 

s+i 



-log^)( S )= £ W^M^ 
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Proof. Note that N(P k ) = N(P) k and hence log N(P k ) = k log N(P) for P G Prim (T) and fceN. 
Hence, from the definition (4.7), we have 

-logZ r (,)- ^ 2^2^^ (logi V(p))m-l 

PePrim(r)n=Ofc=l v b v " 



N{P)~ ks 1 



^ fc m (log iV(P))"^ 1 1 - MP) 
PSPrim (r) k=l x & v " y ' 

logJV(P) iV(P fc )" s+ 5 



P e p£(r)S ^¥ - iV(P fe )^ (log^(^))" 1 ' 

Writing P fc = 7, one obtains the desired expression. □ 

The poly-Selberg zeta function satisfies a differential ladder relation. As a consequence, when 
m > 2, one can obtain an iterated integral representation of Z^ n \s) which gives an analytic 
continuation beyond the line Re (s) = 1. 

Proposition 4.9. (i) It holds that 

(4.9) __l 0gZ M (s) = i og Zr (s) (Re (a) > 1). 

(ii) Fix a G C uni/i Re (a) > 1. T/ien, /or m >2, we have 

f ••• / 2 logZ r (6Ki---^m-i) 

J a J a J a I 

S v ' 7 

m— 1 

where Q ( f l \s,a) := nfe^^r™ *^ («) fcl ( s_a )* . If j $ Spec(Ar), then this gives an analytic 
continuation of Z^ n \s) to the region f2p \ (— 00,— 1], otherwise, to the region Ojt := f2p H {s G 
C|Re(s) > i}. 

Proof. We first show the equation (4.9). The case m = 1 is clear because Zp^(s) = Zr(s). Let 
m > 2. From (4.8) (or the differential equation ^Li m (z) = -Li m -i(z)), we have 

(4.11) ^logZ( m )( S ) = -logZ( m - 1 )( S ). 

Using this equation repeatedly, one obtains the equation (4.9). We next show (ii) by induction on 
m. Assume that | ^ Spec(Ar). Let m = 2. Then, integrating the equation (4.9) with m = 2, we 
have 

(4.12) log zf\s) = log 4%) " f log Z r (0^- 

J a 

Here, we take the path in Re(s) > 1. This immediately shows the equation (4.10) with m = 2 
for Re (s) > 1. Here, in the righthand-side of (4.12), one can move s freely in the region in where 
logZr(s) is single-valued and holomorphic, that is, Or \ ( — 00,— 1] (notice that the Selberg zeta 
function Zr(s) has zeros at s = 1,0, —k for k G N and s = for j G N). Hence the equation 



Milnor-Selberg zeta functions and zeta regularizations 



29 



(2) 

(4.10) gives an analytic continuation of Z-^ ! (s) to f2r \ ( — oo, —1]. Now, assume that the claim (ii) 
holds for m — 1. Then, we have 

m—1 , ^j, 

(4.13) logZ^- 1 )^) = £ ^(e-a^log^- 1 -^) 

fc=0 



(-l) m ^ 2 / / • • • / log Z r (£i)d£i • • • dU-2- 

J a J a J a 



m-2 

Taking the integral of the equation (4.13) together with the formula (4.11), we have 

m— 3 / -. 

'W/„\_1„„7H/„1 i („ „\k+l i_ 7 (m-(fc+l)). 



logZ^( S ) = log 4^ (a) + g ^TJT (* " a) k+1 log4 m - (fe+1 »(a) 

+ (-l) m ~ 2 / / / • • • / logZ r (£iRi • • • dem- 2 de. 



m-2 

Therefore, by the same discussion as above, one shows the claim (ii) for all m > 2. The proof works 
similarly for the case \ E Spec(Ar) (notice that Op - is the connected component of J7r which 
includes the region {s G C | Re (s) > 1}). □ 

Remark 4.10. The discussion above can be applied to a general scheme for arbitrary number fields. 
In fact, let L(s,n) be a L-functions attached to irreducible cuspidal automorphic representation ir 
of GLd(A/f ), where K is an algebraic number field and A^ is the adele ring of K. We have shown 
in [WY] that, as a generalization of the result in [D], a "higher depth regularized products" of the 
non-trivial zeros of L(s, tt) can be evaluated as a product of the Milnor gamma functions and "poly 
L-functions" , which are similarly defined by an Euler product as the poly-Selberg zeta function. 

4.3 Euler product expressions of Z r r (s) 

The following formula can be regarded as an Euler product expression of Zy^(s) (remark that the 

T 

Theorem 4.11. It holds that 



poly-Selberg zeta function Z^ n \s) is defined by the Euler product (4.7)). 



r-i , 1+ v \ (r-m-iy-' 

(r — 1 + m ) ! 



(4.14) Zr.r(«) = " [ z£ +m \s) ^-i-™» {2B - 1) 



\m=0 



Proof. We start from the definition (2.12). Let s <E U + := U n {s G C | Re (s) > |}. Then, using 
the formula 

/ ^ e -K+|)^£ = 2(-Y K w (2ah^) (Re(a) > 0, Re (6) > 0), 
Jo € v o/ 

where K u (x) is the modified Bessel function of the second kind, we have 

e Tr (w,t) 



y/TrT(w + 1 — r) 



x 



76Hyp (r) 
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Here, employing the asymptotic formulas 

1 



(-ip(r-l)!«)|0(4 



T(w + 1 - r) 



and 

K w+ i_ r (tlogN( 1 ))=Ki_ r (tlogN( 1 ))+0(w) 
as w — > together with the equation (see, e.g., [EMOT]) 

'r^i.-i/^n.A-m (r-l + rr 
m=0 



2 2 V2y/ m!(r — 1 — my. 

one can see from (4.8) that 



(4.15) logZ r , r (a) = --^ei>(tM)| 

r— 1 

= (-l) r - x (r - 1)! V ( , r ~ 1 1 +m)! ,, (2tr 1 -" 1 logZ( r+m) ( S ). 
^ mUr - 1 - m) 

m=0 v ; 

This shows the expression (4.14). □ 
Example 4.12. Let t = s - \. Then, we have 

Z T , l (s) = z£\s) = Z r (s), 

z rt2 (s) = z?\s)-WzW( s y 2 , 

Z T>3 (s) = zP{s) 2 ^ 2 Z^{s) 12 ^zf\ S f\ 

Remark 4.13. We remark that, when r > 2 and j £ Spec(Ar), from Proposition 4.9, the 
expression (4.14) also gives an analytic continuation of Zr, r (s) to the region fir \ (— oo, —1]. (We 
have already known this fact from Theorem 4.1. See Figure 3 in Section 1.) 

Finally, as is the case of Zy m (s), we show that the Milnor-Selberg zeta function Zr, r (s) also 
satisfies a differential ladder relation and hence has an iterated integral representation, which again 
gives an analytic continuation to fir \ (—oo, —1] if \ ^ Spec(Ar). 

Corollary 4.14. (i) It holds that 

(4.16) (-^— ^Y' 1 log Z T , r (s) = (r-l)!logZ r (a) (s G U + ). 



(r-l)I 



,2s - 1 ds. 

(ii) Fix a £ U + . Then, for r > 2, we have 



(4.17) Z r , r (s) = Q T , r (s, a) exp (f f' - (jj( 2 " *)) lo § MbHi ■ ■ ■ d£ r -i 



r-1 



w/iere Qv, r (s, a) := Ilfc=o ^r,r-fc(a)^ fel '' (s a ' )k ( s+a 1)fc . // 1 ^ Spec(A r ), i/ien this gives an analytic 
continuation of Zr, r {s) to the region fir \ (— oo,— 1], otherwise, to the region fip~. 
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Proof. To prove the equation (4.16), it is sufficient to show for r > 2 that 

(4.18) (_L_iL)iogZr, r (iO = (r-l)IogZ r , r -i(a). 

Actually, from the expression (4.15) together with (4.11), one sees that 

±logZ r , r (s) = (-l) r (r - 1)1 £ { u + m ~ 2) U 2s ~ ir^log^" 1 ^) 
as ' m r — m — 2)! 

m=0 v ; 

= (2s-l)(r-l)logZ r , r -i(a). 

This shows the equation (4.18). The iterated integral representation (4.17) and the analytic con- 
tinuation can be easily derived by the similar discussion performed in Proposition 4.9. □ 

Acknowledgement. The authors would like to thank the referee and Hirotaka Akatsuka for their 
careful reading of the manuscript and helpful advice. 
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